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Abstract 


The sensitivity of tropical cyclones (TCs) to changes in parameterized convection is 
investigated to improve the simulation of TCs in the North Atlantie. Speeifically, the impact of 
reducing the influence of the Relaxed Arakawa-Schubert (RAS) scheme-based parameterized 
eonvection is explored using the Goddard Earth Observing System versions (GEOS5) model at 
0.25° horizontal resolution. The years 2005 and 2006 characterized by very active and inactive 
hurricane seasons, respeetively, are selected for simulation. 

A reduction in parameterized deep convection results in an increase in TC activity (e.g., TC 
number and longer life cycle) to more realistie levels compared to the baseline control 
configuration. The vertical and horizontal structure of the strongest simulated hurricane shows 
the maximum lower-level (850-950hPa) wind speed greater than 60 m/s and the minimum sea 
level pressure reaching ~940mb, corresponding to a category 4 hurricane - a category never 
achieved by the control configuration. The radius of the maximum wind of ~50km, the location 
of the warm core exceeding 10°C, and the horizontal compactness of the hurricane center are all 
quite realistic without any negatively affecting the atmospherie mean state. 

This study reveals that an increase in the threshold of minimum entrainment suppresses 
parameterized deep conveetion by entraining more dry air into the typieal plume. This leads to 
cooling and drying at the mid- to upper-troposphere, along with the positive latent heat flux and 
moistening in the lower-troposphere. The resulting increase in conditional instability provides an 
environment that is more conducive to TC vortex development and upward moisture flux 
convergenee by dynamically resolved moist eonvection, thereby increasing TC activity. 
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1. Introduction 


This article is inspired by a recent research on tropical cyclone (TC) simulation coordinated 
by the US CLFVAR hurrieane working group r http://www.usclivar.org/working-groups/hurricane') 
Among various science issues raised in the research, it was found that many current general 
eireulation models (GCMs) seriously underestimate tropieal cyelone (TC) aetivity over the North 
Atlantic when run at -0.5 degree latitude/longitude or coarser horizontal resolution as opposed to 
some other basins such as the North Pacific (Walsh et al. 2013; Shaevitz et al. 2014). Only a few 
recent GCMs simulate TC activity in the North Atlantic with some success (e.g., Zhao et al. 
2009). A question is raised on whether horizontal resolution finer than -0.5 degree is necessary 
to achieve reasonable TC numbers over the North Atlantic, even though resolution appears less 
critieal for simulating TC numbers than the intensity (Straehan et al. 2013). Compared to TC 
numbers, simulating realistic TC intensity and spatial structure appears more unattainable at 
these resolutions (-0.5 degree or eoarser), with generally weak and poorly organized systems 
occurring in all the GCMs (more than 10) involved in the working group. Higher resolution may 
be more appropriate to allow for a better representation of Afriean easterly waves and the 
associated strong low-level vorticity, leading to more realistic simulation of Atlantic TC activity 
including number, intensity, life-cycle and spatial structure (Knutson et al. 2007; Caron et al. 
2011; Putman and Suarez 2011; Daloz et al. 2012; Emanuel 2013). 

Several previous studies have identified deficieneies in the simulation of moist convective 
processes as substantially affecting the ability to simulate TC characteristics (Slingo et al. 1994; 
Zhang and MeFarlane 1995; Smith 2000; Benjamin et al. 2008; Zhao et al. 2012). Smith (2000) 
suggests that major improvements in TC predictions will depend on improvements in the 
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representation of convection in hurricane models. Knutson et al. (2004), Atlas et al. (2005), 
LaRow et al. (2008), Ma and Tan (2009) and Reed and Jablonowski (2011) performed model 
experiments to examine the sensitivity of the global TC numbers to changes in the convective 
parameterizations. Shen et al. (2006) and Stan (2012) examined the impact of disabling the 
convective parameterizations (in particular the schemes developed by Arakawa and Schubert 
(1974) and Moorthi and Suarez (1992)) on GCM hurricane forecasts. Both studies agreed on that 
explicit representation of cloud processes produces a larger number of TC events, with stronger 
intensity and longer life-cycles (Reed and Jablonowski 2011; Stan 2012). However, the details of 
the atmospheric processes responsible for altering TC activity were not the focus of the above 
studies. Some of the atmospheric responses to changes in deep convective activity are discussed 
in Zhao et al. (2012), which focuses on the impacts on vertical velocity, temperature and 
humidity. That study identifies stronger sensitivity of mid-tropospheric vertical motion to the 
minimum entrainment compared with the other variables (humidity and vertical wind shear): the 
importance of the response of vertical motion is also discussed in Oouchi et al. (2006). 

The sensitivity of the atmospheric response to the convective parameterization and its 
impact on TC activity motivates the current study to improve the simulation of TC activity over 
the North Atlantic in GCMs. As in Tokioka et al. (1988), Held et al. (2007) and Kang et al. 
(2008), we find that increasing the threshold of mi nimum entrainment rate makes the convective 
plume entrain more dry air as it ascends, losing in-cloud buoyancy more rapidly when entraining. 
As a result, deep convection is inhibited from occurring in updrafts with a lateral entrainment 
rate 2 lower than a threshold value 2o (Tokioka et al. 1988; Held et al. 2007). Parameterized 
deep convection, including the sub-grid scale cloud representation is restricted by increasing 2o 
(Held et al. 2007; Kang et al. 2008), while the fraction of large-scale condensation by non- 
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convective clouds increases through the large-scale cloud/condensation module (Tiedtke 1993) 
as opposed to the convection module (Kang et al. 2008). Wang (2014) found that this non-deep 
convection plays a role in moistening the lower to middle troposphere whereas deep convection 
moistens the upper- troposphere. Zhao et al. (2012) also found that an inhibition of the convective 
parameterization through enhanced lateral mixing into convective plumes leads to increased TC 
genesis, along with a eolder and drier upper troposphere. 

The main objective of this study is to improve the simulation of North Atlantic TC activity 
in the National Aeronautics and Space Administration (NASA) Goddard Earth Observing 
System Version 5 (GEOS-5) model (Rienecker et al. 2008; Molod et al. 2012), including the TC 
number, intensity, life eycle and horizontal/vertical structures. We focus on the cumulus 
entrainment rate in the GEOS-5 model run at a relatively high horizontal resolution (0.25 degree) 
to determine the thresholds of minimum entrainment that reliably reproduee the TC activity over 
the North Atlantic. An important aspect of this study is that we also focus on maintaining a 
realistic mean elimate state - an aspect of the simulations (the response of the mean state to the 
changes in deep convection) many previous studies have not emphasized. The analysis of the 
runs is geared to better understanding the atmospherie response determining the TC activity over 
the North Atlantic with a particular focus on the atmospheric instability over the TC genesis 
region, the thermodynamie and radiative balance, the low-level fluxes and circulation, and the 
mean climate state response. 

The organization of the paper is as follows. Section 2 describes the GEOS-5 model and the 
experimental design. The results of the sensitivity experiments are presented in Section 3, 
including TC numbers, intensity, life cycle and spatial structure of the strongest hurrieanes. We 
investigate in Section 4 the relevant atmospheric dynamics and physics to explain the TC activity 
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changes that occur as a result of the modifications to the parameterized deep convection. Some 
verifieation of the atmospheric basic state changes due to the deep convection changes is also 
provided in Section 4, followed by concluding remarks and discussion in Section 5. 

2. Model and Experimental design 
2.1 Model 

We utilize the NASA GEOS-5 model (Rienecker et al. 2008; Molod et al. 2012) for our 
experiments. The model is run with 72 hybrid-sigma vertical levels, extending to 0.01 hPa, and 
~0.25° latitude/longitude horizontal resolution. The convection scheme in the GEOS5 model is 
basieally a modified version of the Relaxed Arakawa Schubert (RAS) scheme of Moorthi and 
Suarez (1992). In the modified version, multiple convective plumes in RAS have the convection 
base level at the top of the planetary boundary layer. The adjustment time scale of plumes is a 
function of the vertical depth of plumes, which varies from 0.5 hr for the shallowest plume to 12 
hr for the deepest. Another major modifieation to RAS, which affects convective variability of 
the model substantially and hence the intensity of simulated TCs, is the stochastic determination 
of the minimum entrainment threshold. This is basically same as in Tokioka et al. (1988), which 
modification is known as limiting frequent adjustment of deep convection in a model and 
inereasing the variability of deep convection (Eee et al. 2003; Ein et al. 2008). Eor “Tokioka” 
constraint in this study, the threshold value is non-deterministic so that the subsequent value is 
determined randomly with a predetermined power function: 

x", 0<x<l (1) 

a D 

where a is a free parameter, D the depth of the subcloud layer, and x a random number chosen at 
every model integration time step. Here the depth of the subcloud layer D is assumed to be 
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related with the plume radius by Simpson and Wiggert (1969), and this stochastic determination 
has a basis to consider spatial statistics of the observed cloud widths (or radius) and volumes that 
follow a power law distribution found by Bacmeister and Stephens (2011) from the CloudSat 
data (Stephens et al. 2002). Figure 1 compares the mi nimum entrainment threshold values 
between the conventional, deterministic method from Tokioka et al. (1988) and the stochastic 
method in this study. In case of n=l/2, which is used in this study, the threshold curve increases 
with the square root function (Fig. 1), and restores to the conventional Tokioka constraint. This 
tends to impose restriction on the entrainment rate in most cases, but occasionally allowing a 
very deep and less entraining convection plumes. Putman and Suarez (2011) used the same 
version of GEOS-5 that used in this study for their tropical cyclone simulation, but with much 
higher horizontal resolution of ~7km. 

Increase in the Tokioka constraint (i.e., increase in the threshold of minimum entrainment) 
makes more difficult condition for sub-grid scale deep convection. While the parameterized deep 
convective process is suppressed as a result of this Tokioka constraint, resolved-scale non- 
convective clouds are favored and thus, more active large-scale cloud/condensation and 
precipitation process becomes feasible. 

2.2 Experimental design 

Three types of experiments are carried out by applying different Tokioka constraints in 
model simulations of the years 2005 and 2006, which are very active and inactive hurricane 
years, respectively. We call these three experiments A, B, and C. Experiment A (ExpA) has the 
strongest Tokioka constraint, where the greatest possibility of large thresholds of minimum 
entrainment exists, while Experiment B corresponds to an intermediate range of thresholds 
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(ExpB), and Experiment C (ExpC), which is a baseline control run, corresponds to the lowest 
minimum entrainment range that allows more parameterized deep convection relative to the 
ExpA and ExpB. These were implemented in the model with the minimum a/D = 2 km'^ (ExpA), 
0.7 km'^(ExpB), and 0.45 km'^ (ExpC), respectively in Equation (1). Please note that solid curve 
as a bottom line for minimum entrainments in Eig. 1 can vary with the a/D values given in each 
experiment. The curve for ExpA is placed above the other two curves (for ExpB and ExpC), 
while the opposite is true for the curve for ExpC. Although the choice of n in Equation (1) can 
affect the results, but this is left for future sensitivity study. As mentioned earlier, the horizontal 
resolution is a quarter-degree latitude/longitude. The minimum entrainment thresholds is 
increased in our experiments (e.g., ExpA and ExpB) to identify the model sensitivity in 
simulating TCs with a focus on better performance across the Atlantic basin. Erom the above 
experiments, we intend to estimate the best Tokioka constraint as determined by the 
improvements in TC simulation including genesis, lifetime, track and intensity at this resolution. 
Each type of experiments consists of three ensemble members. All simulations are done with 
prescribed weekly sea surface temperature (SST) forcing (the HadlSST of Rayner et al. 2003), 
and are initialized from Modem-Era Retrospective Analysis for Research and Applications 
(MERRA, Rienecker et al. 2011) data on May 1 each year. 

3. Sensitivity of TC characteristics to changes in parameterized deep convection 

TCs are detected for each year (2005 and 2006) and each experiment (ExpA, ExpB and 
ExpC). The number of detected tropical storms and hurricanes is shown in Table 1. The basic 
algorithm of TC detection and tracking is the same as that used in Vitart et al. (2003), Knutson et 
al. (2007) and EaRow et al. (2008). Threshold values for wind speed, relative vorticity, warm 
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core, sea level pressure and minimum duration time are assigned in the algorithm to define TC 
(Table 2) at a quarter degree resolution based on suggestions in Walsh et al. (2007). The 
threshold values used here can be considered realistic because they are very close to the 
thresholds to detect TCs from observational data. Exp A produces an excessive number of 
tropical storms and hurricanes, compared with observations shown on the rightmost column 
(Table 1). In partieular, it produces too strong TC aetivity (hurricanes more than 20) even in 
2006, which was recorded as a relatively calm hurricane year (Knapp et al. 2010). ExpB 
produces more reasonable TC numbers, as they are closer to observations than either Exp A or 
ExpC. ExpC shows generally weak TC activity and produces much weaker than observed TC 
activity during the aetive hurricane season of 2005. 

TC tracks are plotted for each experiment. During the 2005 season (Eig. 2), a variety of 
traeks were observed; 1) Cap Verde-type TC systems that crossed the Atlantic making landfall 
over the American continent, 2) early recurvers, 3) Gulf of Mexico TCs, and also 4) eastward 
traeking systems embedded in the westerlies. Comparison between the experiments inideates that 
ExpA and ExpB both display a variety of tracks much larger than ExpC, and closer to the 
observation. The results for 2006 (Eig. 3) show features similar to those obtained for 2005 in the 
sense that the increased thresholds of minimum entrainment produce larger number of TCs with 
longer life cycles and a track variety more similar to the observation over the Atlantie. This 
model response is consistent with Emanuel et al. (2008) and Stan (2012), that showed a 
sensitivity of TC aetivity to changes in the cumulus parameterization in the Community Climate 
System Model (CCSM) (Gent et al. 2011). 

Additional details of the model responses (including TC activity, TC intensity and vertical 
hurricane structure) produced in each experiment are shown in Eig. 4 through 8. Scatter plots of 
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low-level (850-950hPa) wind versus SLP for 2005 (Fig. 4), show that the maximum wind speed 
is up to 61ms'' in the ExpA, 60ms'' in ExpB, and less than 60ms'' in ExpC. The minimum SEP 
is down to 936, 943, and 963hPa, respectively in the three experiments. ExpA and ExpB reveal 
improved distributions of SEP and wind speed compared with ExpC. Note however that the 
observed wind speed is up to 80ms'' and minimum SLP deeper than 920hPa for extremely strong 
hurricanes - something not reproduced by the model in any of the experiments. Much finer 
horizontal resolution (finer than 10km) such as that of Noda et al. (2012) and Putman and Suarez 
(2011) appear to be necessary for capturing such extremely strong hurricanes. The maximum 
wind speeds identified in ExpA and ExpB are comparable to those found in the idealized 
hurricane experiment of Reed and Jablonowski (2011). They simulated maximum wind speeds 
of ~60ms'' at ~lkm height and a quarter-degree horizontal resolution. 

Eigure 5 is the same as Eig. 4 but for 2006. The results show that the simulated wind speeds 
and minimum center SLP deepening for very strong hurricanes is comparable to or stronger than 
those for the observations (ExpA and ExpB). In contrast, ExpC exhibits weaker maximum wind 
and minimum SLP distribution compared with the observations. 

The vertical structure of the strongest hurricane for ExpA is plotted in Eig. 6. The top left 
panel shows the wind speed (shaded) and temperature (contoured), and the right panel shows the 
warm core. The results for 2005 show a compact core and hurricane eye that is well defined. The 
maximum low level wind speeds (greater than 60ms'') and the radius of maximum winds (less 
than 50km) both correspond to very realistic TC representation. The warm core value is greater 
than 10°C and situated in the upper-troposphere, quite similar to the typical structure of observed 
hurricanes (Erank 1977). The bottom panel is the time series of the minimum SLP and maximum 
wind speed following the moving hurricane. The minimum SLP drops to 935hPa and maximum 
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wind speed reaches 61ms'\ The main features for 2006 are quite similar to those found for 2005, 
showing a well-defined hurricane characterized by a maximum wind speed of 65ms'' and 
minimum SLP of 935hPa. 

The vertical structures of the strongest hurricanes in ExpB are shown in Fig. 7. Compared 
with ExpA, they are a little weaker, but still show a reasonable vertical structure in terms of 
compactness of the core, with the radius of the maximum wind less than ~50km (2 grid points in 
longitudinal direction), the maximum wind value near 60ms'', the minimum SEP of about 
940hPa, and vertical profile of the warm core (reaching 10°C), which corresponds to a category 
4 hurricane based on the Saffir-Simpson scale ( http://www.nhc.noaa.gov/aboutsshws.php) in 
2005. Observed wind .vs. SLP scatter plots in Figs. 4 and 5 indicated that the strongest observed 
hurricane in 2006 is weaker than that in 2005. The model reproduces well this difference 
between the two years (Fig. 7), but the difference is not as distinguishable as that found in the 
observation. As discussed in Fig. 4, the observed category 5 hurricane in 2005 was very hard to 
reproduce through a quarter-degree resolution run, resulting in rather a small difference in the 
intensity of the strongest hurricane between the two years in Fig. 7. 

ExpC, which is the control run, also shows reasonable results (Fig. 8), though the storms are 
less organized than those found in the other two experiments. The warm core is weak with a 
value of approximately 6~8°C and located in the mid-troposphere. The wind speed, SLP, and 
sharpness of the hurricane core are also not as realistic as those in ExpA and ExpB. The weak 
hurricane structures found in ExpC are not inconsistent with the study of Vitart et al. (2001), 
which found that the vertical structures of hurricanes simulated by RAS tend to be less intense 
and have lower warm cores than those produced with moist convective adjustment (MCA) 
schemes (Manabe 1969), though the RAS scheme produces more TCs and higher relative 
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humidity than MCA schemes. In any event, the hurricane vertical structures shown in Fig. 6 
through 8 provides eompelling evidenee that the suppression of parameterized deep eonveetion 
tends to produee more intense hurrieanes with better organized vertieal structures. 

4. Atmospherie response to parameterized deep eonveetion 

The physieal meehanisms by whieh the suppression of the eumulus parameterization 
promotes TC activity are investigated next. The role of atmospheric moist static stability - a 
factor kn own to affect TC characteristics (Smith 2000) is shown in Fig. 9, where the vertieal 
profile of moist static energy (left panel) and the vertical gradient of equivalent potential 
temperature (right panel) are eomputed over the TC genesis region, whieh eovers most of the 
eentral and eastern tropieal Atlantic (60°W~15°W, 5°N~20°N). Profiles from ExpA, ExpB, and 
ExpC all exhibit a deerease with height in moist statie energy at low levels (below 700~800hPa 
height) and an inerease at upper-levels, implying the possibility of eonditional instability. The 
ExpA result shows the possibility of stronger instability than ExpB and ExpC. The vertieal 

gradient of equivalent potential temperature (— (Eig. 9b) also shows that lower tropospherie 

vertieal column with the negative is more extensive in Exp A than the other two 

experiments, indicating the most unstable atmosphere occurs in ExpA, followed by ExpB and 
ExpC. 

The vertical temperature and humidity profiles are plotted over the TC genesis region in Eig. 
10 to indicate how they are linked to the atmospherie instability. Here we plot the differenees in 
the vertical profiles of temperature/humidity between two experiments. The most striking 
feature of these results is the mid- to upper-tropospherie cooling (Eig. 10a) and drying (Eig. 10b) 
and lower- tropospheric moistening (Eig. 10b) resulting from the reduction in parameterized deep 
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convection. This vertical structure of the humidity is in good agreement with the results of Zhao 
et al. (2012). With a larger threshold of minimum entrainment, it is not unexpected that more dry 
air will be entrained into the plume, resulting in a drier upper atmosphere with decreased 
buoyaney (Held et al. 2007). The reduction of in-cloud buoyancy leads to the suppression of 
parameterized deep convection and less upper-tropospheric condensation heating. Figure 10 
reflects such a change in the vertieal profile of temperature and moisture, providing favorable 
conditions for developing an unstable atmosphere. Upward moisture flux and dynamically 
resolved convection are more likely to occur after the lower troposphere is moistened enough 
(Wang 2014) under this unstable atmospheric condition. 

Figure 1 1 shows the latent heating and evaporative flux with a foeus on the near surface 
moisture flux. We note that the differences between ExpB and ExpC exhibit very similar patterns 
though with smaller magnitude (Eigure not shown). Eigure 11a shows that the latent heat flux is 
significantly increased in ExpA throughout the Atlantic basin. The enhanced evaporative flux 
from the ocean (Eig. 1 lb) gives rise to the increase in latent heat flux. We will see later (Eig. 12) 
that the upward flux at the surface is also enhanced by increased ascending motion in the 
unstable atmosphere. This is consistent with Zhao et al. (2012), who found that surfaee heat 
fluxes including latent heat flux were influenced by vertical motion, and that this was important 
for making the resolved-scale convection more favorable for tropieal cyelone genesis. 

Eigures 11c quantifies the increase in conditional instability we inferred from Eig. 9 through 
11b, in terms of the convective available potential energy (CAPE). It shows that CAPE is 
significantly enhanced in ExpA (compared with the control) through the constraint on the 
parameterized deep convection. CAPE exists within the conditionally unstable layer of the 
troposphere, when a lifted parcel of air is warmer than the ambient air. It can be interpreted. 
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based on Fig. 10, that strong upper-tropospheric cooling overlying relatively warmer air at 
lower-levels accompanied sufficient moisture, contributes to a build up of CAPE. The larger 
CAPE in ExpA compared with the control run (ExpC) also demonstrates that dynamically 
resolving convective processes would be enhanced with the suppression of parameterized deep 
convection, resulting in an increase in large-scale condensation and precipitation (see the 
fractional precipitation change in Eig. 15). 

Eigure 12a clearly illustrates the strong response of the vertical motion (see also Oouchi et al. 
2006) to the unstable atmosphere caused by the restriction of parameterized deep convection. 
The increase in atmospheric instability over the TC genesis region leads to increased ascending 
motion throughout the vertical column with the largest increase occurring in ExpA, followed by 
ExpB and then ExpC. The enhanced vertical ascent is associated with enhanced low-level 
convergence (Eig. 13). This low-level convergence, combined with low-level moistening (Eig. 
10) and positive latent heat flux (Eig. 11a), plays a role in increasing the moisture flux 
convergence at low levels (Eig. 12b), with the largest values occurring in ExpA, followed by 
ExpB and ExpC. 

Eigure 13 shows the distribution of the differences between ExpA and ExpC in some other 
atmospheric quantities. The results show that low-level vorticity (850hPa), near-surface wind 
variance (10m), low-level convergence (925hPa) and vertically integrated moisture flux 
convergence are all enhanced over the eastern Atlantic. The enhancement is more pronounced 
over the TC genesis region between about 10°N and 15°N (for genesis latitudes, see also Eigs. 2 
and 3). These spatial patterns emphasize once again that, as the parameterized deep convection is 
suppressed through an increase in the minimum entrainment threshold, the resulting increase in 
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atmospheric instability makes dynamically resolved moist convection more active by enhancing 
low-level atmospheric motions (e.g., convergence, vorticity and wind variance). 

To quantify that the atmospheric structures discussed so far provide a favorable TC genesis 
environment, we next calculate the genesis potential index (GPI). Here, the original version of 
GPI (Emanual and Nolan 2004) has been modified following Murakami et al. (2011) to include 
the influence of vertical motion in contributing to TC genesis. The modified GPl' is given by 

GPI=|10=,|«^ (1 + o.iv;)-^ (^), (2) 

where rj is the absolute vorticity at 850hPa level, RH is the relative humidity in percent at 

850hPa, Vpot is the maximum potential intensity (MPI; m s'^) defined in Emanuel (1995), Vg is 
the vertical wind shear (E(850hPa) minus E(200hPa)), and w is the vertical wind velocity (Pa s'^) 
at SOOhPa. Epot is modified by Bister and Emanuel (1998) from the original version of Emanuel 
(1995) so that 

^ - CAPE^), (3) 

where C* and Cd are the exchange coefficient for enthalpy and the drag coefficient, respectively. 
Ts is SST and To is the mean outflow temperature. Two CAPE-related quantities, CAPE and 
CAPE^ are the CAPE of the air displaced upward from saturation at sea level with reference to 
ambient air and the CAPE of the air at boundary layer, respectively. 

The GPI formulation considers the influence of all the horizontal and vertical atmospheric 
structures and atmospheric instability discussed in reference to Pig. 9 through 13. As shown in 
Pig. 14, the GPI exhibits larger values in ExpA and ExpB than in the control run (ExpC). Earge 
' The source code for the PORTRAN version of the GPI calculation is available from 
ftp://texmex.mit.edu/pub/emanuel/TCMAX/pcmin revised.f 
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values occur over the Atlantic in low-latitudes, the Caribbean Sea, and Gulf of Mexico, where a 
majority of TC genesis takes place. 

An important consideration is whether the above modifications in the convective 
parameterization produce changes in the atmospheric mean state and its variability (e.g., seasonal 
cycle) that are within acceptable ranges. Slingo et al. (1994) and Zhang and McFarlane (1995) 
found some sensitivity of the tropical mean climate to the parameterized moist convection. 
Figure 15 shows the differences between each experiment in the mean precipitation for 2005. We 
plot both the large-scale precipitation and convective precipitation to indentify the fractional 
precipitation changes. The left panels show the difference between ExpA and ExpC in total 
precipitation (Eig. 15a), large-scale precipitation (Eig. 15b) and convective precipitation (Eig. 
15c). The results show that both total and large-scale precipitation increase noticeably over the 
tropics for ExpA, whereas there is a substantial reduction in convective precipitation. In contrast, 
for ExpB (the right panels) there is only a very small change in total precipitation over the 
tropics (Eig. 15d), while the large-scale precipitation over the TC genesis region increases 
slightly (Eig. 15e) along with a decrease in convective precipitation (Eig. 15f). This fractional 
precipitation change is due to an increased activity of grid-scale cloud/condensation (i.e., large- 
scale precipitation) along with a suppression of parameterized deep convection (i.e., convective 
precipitation) (Reed and Jablonowski 2011; Zhao et al. 2012). 

The area-averaged difference in total precipitation over the Atlantic basin between ExpB 
and ExpC (Table 3) is 0.03 mm day’' in 2005 and 0.09 mm day’' in 2006, while for ExpA the 
precipitation increases with respect to ExpC by 0.5 and 0.6 mm day'^ in those two years, 
respectively. All three experiments overestimate the observed precipitation (left column in Table 
3; see also Eig. 16a). Since the observationally-based precipitation estimates from the Global 
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Precipitation Climatology Project (GPCP) version 2.2 (Adler et al. 2003; Huffman et al. 2009) 
are provided on a relatively coarse grid (2.5° lon.-lat.), we speculate that part of the apparent 
overestimation by the model might in fact represent an underestimation of the coarser resolution 
GPCP precipitation. Table 3 indicates that the excessive total precipitation in ExpA makes it a 
less than ideal model to use even though it has more intense hurricanes. Also, ExpA produces an 
excessive number of storms during the inactive hurricane year of 2006 (Table 1). In contrast, 
ExpB seems a better option since it has minimal impact on the total precipitation (amount and 
spatial distribution), and it has improved 3-dimensional hurricane structures, intensities, and 
numbers, compared with the control (ExpC). 

Eigure 16 compares for each experiment the seasonal variation (June through November) of 
precipitation, SEP, surface temperature and upper-level wind speed over the typical TC genesis 
region. We see that ExpB tends to remain close to ExpC in all the variables, while ExpA shows 
a noticeable increase in total precipitation and surface temperature, and decrease in SEP 
compared to the control run. A comparison with observations indicates that, except for the upper- 
level wind speed, ExpA tends to have larger biases than either ExpB or the control (ExpC). 
Overall, our results suggests that a small reduction in the influence of the RAS convective 
parameterization (ExpB) produces improved North Atlantic TC activity without having a 
negative impact on the mean climate of the model when run at Vi degree spatial resolution. 

5. Concluding remarks and discussion 

This study investigated the sensitivity of the North Atlantic TC activity to changes in 
parameterized deep convection in the NASA GEOS-5 model run at 0.25 degree 
latitude/longitude horizontal resolution. The study found that reduction in the influence of the 
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RAS scheme resulting from an inerease in minimum entrainment for parameterized deep 
eonvection improves/enhances (eompared to the standard eontrol model) the TC aetivity in terms 
of numbers, intensity, life eycle, and the 3 -dimensional storm struetures. Our case studies for the 
2005 (very aetive) and 2006 (very inaetive) hurrieane seasons revealed that only a modest 
inerease in minimum entrainment (ExpB) was neeessary at a quarter-degree resolution to 
produee TC numbers reasonably elose to observations. The strongest hurrieane reaehed a 
minimum SLP of 940hPa and low-level maximum wind speeds greater than 60m/s. The vertieal 
strueture is eharaeterized by a well-defined hurrieane eye, an upper-tropospherie warm core 
value exeeeding 10°C, and maximum winds loeated at low levels within 50km of hurricane 
eenter. These magnitudes are comparable to, or improvement over, those of Shen et al. (2006) 
and Reed and Jablonowski (2011) at similar resolution. Furthermore, we found that a stronger 
eonstraint on the parameterized deep eonveetion (i.e., ExpA) produeed even stronger hurrieanes 
(with minimum SEP of 935hPa, and a maximum wind speed of 65m/s), but this strong reduetion 
in parameterized eonveetion was not desirable as it produeed strong TC aetivity even in a weak 
hurricane year (2006) and lead to inereased bias in the atmospheric mean state. 

Our analysis of the causes of the TC ehanges suggests that the following proeesses are at 
work. First, an inerease in the threshold of minimum entrainment eauses more entrainment of dry 
air into the eonveetive plume. This entraining proeess leads to a reduetion in buoyaney, resulting 
in less favorable conditions for deep eonvection. This leads to a reduction in condensation 
heating in the upper-troposphere, resulting in upper-tropospherie eooling. As a result, the upper- 
troposphere becomes drier and cooler, while the lower-troposphere beeomes moister and warmer. 
This is aeeompanied by an inerease in positive near-surfaee latent heat flux. This ehange in the 
vertieal structure of the atmosphere is apparent over the Atlantic basin, easing the condition for 
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developing atmospheric instability. The greater instability is evidenced by changes in the profile 
of moist static energy, the vertical gradient of equivalent potential temperature, and the 
distribution of CAPE. As a consequence to the changes in atmospheric instability there is an 
increase in a) low level positive vorticity and convergence, b) upward motion, and 3) low-level 
moisture flux convergence. These changes help the development of explicit-scale convection and 
subsequent large-scale condensation of wave vapor, and represent a transition of the moist deep 
convection from the parameterized kind to that occurring at the resolved scales. This is reflected 
by a decrease in convective precipitation contributed by parameterized convection and an 
increase in large-scale precipitation. The above interpretation is supported by the changes in the 
GPI (Emanuel and Nolan 2004; Murakami et al. 2011) showing (for ExpB) an enhancement in 
the TC activity over the Atlantic, with more and stronger TCs than the control run (e.g., 
category 4 hurricanes were produced, whereas the control never produced hurricanes stronger 
than category 3). 

The explicit-scale moist convection favored by atmospheric instability appears to 
compensate for the upper-tropospheric cooling and drying produced by the reduction in 
parameterized deep convection, with the atmosphere acting to maintain thermal equilibrium in 
the upper-troposphere. This likely explains why, at least for the experiment with only a modest 
change in the threshold of minimum entrainment (ExpB), the total precipitation remains close to 
those of the control. 

Due to the limited number of cases investigated here, it is premature to say that ExpB truly 
represents an optimal setting for the reliable simulation of TCs at this resolution. Nevertheless, 
the current set of experiments are encouraging in that they demonstrate how years with widely 
differing TC activity (2005 and 2006) can be reproduced by global climate model with a rather 
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minor change to the convective parameterization, and that this can be done without degrading the 
mean state. 

This article suggests that quarter-degree horizontal resolution may be the minimum 
resolution necessary for achieving realistic TC simulations via modifications to the 
parameterized convective process (Reed and Jablonowski 2011), while still maintaining a 
realistic mean state. This eonelusion is based in part on other experiments we earried out at half 
and one-third degree resolution (not shown), in which realistic TC numbers were difficult to 
achieve without modifying the threshold values for TC deteetion. Of course the results at Vi 
degree resolution still have several limitations. Key among them is the inability to produce the 
most intense (category 5) hurrieanes. It appears that even at substantially finer than a quarter- 
degree resolution, fine scale hurricane structure can be improved through the modification of the 
convective parameterization. For example, Putman and Suarez (2011) concluded from the 
analysis of GEOS-5 forecasts run at 7-km horizontal resolution, that intense precipitation was 
better formed within the hurricane eye-wall and surrounding rain bands with a reduetion in the 
influence of the RAS convective parameterization, compared to a control in which RAS was not 
modifed. As such it appears that, even at considerably higher resolution than considered here, 
improvements in the simulation of TC activity can still benefit from a proper modification of the 
convective parameterization. 
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588 Table 1 The number of TCs (middle and bottom row) detected for each year and experiment. 

589 Rightmost column represents the observed TC numbers archived at International Best Track 

590 Archive for Climate Stewardship (IBTrACS) (Knapp et al. 2010). The values in parenthesis in 

591 the middle and bottom row are the number of hurricanes. 



ExpA 

ExpB 

ExpC 

(control run) 

Observation 

2005 

36 (24) 

25 (15) 

18(8) 

29(15) 

2006 

27 (21) 

13(8) 

8(6) 

10(5) 
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594 

595 Table 2. List of the threshold values for deteeting the TC using the deteetion/tracking algorithm 

596 based on Vitart et al. (2003) 


Variables 

loeal relative 
vortieity 
maximum 
(850hPa) 

warm core 
temperature 
maximum 

minimum sea 
level pressure 
(SLP) 

minimum 
lower-level 
wind speed 

minimum 

duration 

Criteria 

1.0 xlO'^s'' 

6K, 

Distance 
between the 
TC center and 
the center of 
the warm core 
must not 
exceed 2° 
lon.&lat. 

Minimum 
SLP defines 
the TC center 
and must 
exist within 
2°x2° radius 
of the 
vortieity 
maximum 

17m s'^ 

4 days 
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602 


603 Table 3. Difference in total precipitation averaged over the North Atlantic basin that covers 

604 90°W-0°E and 0°N-50°N. Difference between ExpA and control run is shown in the middle 

605 column whereas the difference between ExpB and control run in the right column. The second 

606 and the third row represent the result for the year 2005 and 2006, respectively. Observed 

607 precipitation on the left column is obtained from the Global Precipitation Climatology Project 

608 (GPCP) dataset, which is 2.5 lon.-lat. gridded. 



ExpA minus ExpC 

ExpB minus ExpC 

2005 (Obs.=3.62) 

4.89-4.40=0.49 

4.43-4.40=0.03 

2006 (Obs.=3.57) 

5.02-4.42=0.60 

4.51-4.42=0.09 
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612 

613 

614 Figure 1 Profile of Xq as a funetion of x in equation (1) with n=l/2 (solid eurve). A straight 

615 dashed line represents the deterministie A,o=a/D, which applies the Tokioka constraint identically 

616 at every time step. The entrainment rate of convective plumes need to be higher than ko, which 

617 range is represented with the black hachured area for the deterministic method. The entrainment 

618 rate allowed for the stochastic method is represented with both the black and grey hachured areas. 

619 
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Atlantic TSC tracks (2005) 



620 

621 

622 Figure 2. TC tracks for 2005 for a) ExpA, b) ExpB, c) ExpC and d) Observations. 
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Atlantic TSC tracka (2006) 
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Figure 3. Same as Figure 2, but for 2006. 
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Figure 4. Scatter distribution of sea level pressure (SLP) and lower-level wind speed of TCs in 
2005 for a) Exp A, b) ExpB and c) ExpC. 
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SLP <& Wind scatter diagram (2006) 
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637 Figure 5. Same as Figure 4 but for 2006. 
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Figure 6. a) and b): Longitude-height eross-section of the strongest hurrieane in 2005 simulated 
by ExpA. Shaded and contoured in a) is wind speed and temperature, respectively, and b) 
represents warm core computed as temperature deviation from zonal mean over 40° longitudes 
with a hurricane core centered. Black line on the bottom panel denotes the time evolution of the 
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maximum wind speed at the lower-level (850-950hPa) of hurricane center whereas blue line the 
minimum SLP. 
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Figure 6: (Continued) Same as the figure on the previous page but for 2006. 
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Figure 7. Same as Figure 6 but for ExpB 
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Figure 7: (Continued) Same as the figure on the previous page but for 2006. 
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Figure 8. Same as Figure 6 but for ExpC. 
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Figure 8: (Continued) Same as the figure on the previous page but for 2006. 
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Figure 9. Vertical profile of area-averaged a) moist static energy (J/Kg) and b) vertical gradient 

dQ 

of equivalent potential temperature from ExpA (black), ExpB (blue) and ExpC (red) 

experiments for 2005. Geographical domain for area-averaging covers the typical TC genesis 
region over the Atlantic that spans 60°W-15°W, 5°N-20°N. 
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Figure 10. Left: Vertical profile of area-averaged difference in a) temperature and b) specific 
humidity for 2005. Black curve denotes difference between ExpA and ExpC (ExpA minus 
ExpC) while the blue curve denotes difference between ExpB and ExpC (ExpB minus ExpC). 
Geographical domain for area-averaging covers the typical TC genesis region over the Atlantic 
that spans 60°W-15°W, 5°N-20°N. 
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Figure 11. Difference in a) surface latent heat flux, b) evaporative flux, and c) convective 
available potential energy between ExpA and ExpC (ExpA minus ExpC) during hurricane season 
(June through November) of 2005. 
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Figure 12. The vertical profile of area-averaged a) omega velocity and b) moisture flux 
convergence for the hurricane season (June through November) of 2005. The geographical 
domain for area-averaging covers the typical TC genesis region over the Atlantic that spans 
60°W-15°W, 5°N-20°N. Black, blue and red curves denote results for ExpA, ExpB and ExpC, 
respectively. 
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Figure 13. Difference in a) relative vorticity at 850hP, b) wind variance at 10m level, c) 
convergence at 925hPa and d) vertically integrated moisture flux convergence between ExpA 
and ExpC (ExpA minus ExpC) during hurricane season (June through November) of 2005. 
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Figure 14. Distribution of modified genesis potential indiees for 2005 from a) ExpA, b) ExpC 
and e) difference (ExpA minus ExpC) 
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Figure 15. Difference in a) total precipitation, b) large-scale precipitation and c) convective 
precipitation between ExpA and ExpC (ExpA minus ExpC) during 2005 hurricane season. Right 
panel is the same as the left panel but for difference between ExpB and ExpC (ExpB minus 
ExpC). 
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Figure 16. The seasonal variation of the area-averaged a) total preeipitation, b) SLP, c) 2m air 
temperature and d) 200hPa wind speed for 2005. The geographical domain for area-averaging 
covers the typical TC genesis region over the Atlantic that spans 60°W-15°W, 5°N-20°N. Black, 
blue, and red solid lines represent the results from ExpA, ExpB and ExpC, respectively. The 
observed (GPCP) precipitation and remaining variables from MERRA are plotted with black 
dashed lines. 
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